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Tij- ' Abstract 

(N : „ 

We calculate the full O(a^a) corrections to the process of gluino pair production 
at hadron colliders in the framework of the real MSSM. We show that these 
Qh' contributions can be neglected at the LHC performing a scan over a wide region 

of the parameter space. The impact of these corrections in the parameter range 
O I investigated at the Tevatron is small. 

> ; 

2° ; 1 Introduction 

m ■ 

I Despite of its success, the Standard Model (SM) of particle physics is affected by 

OO ■ theoretical and phenomenological problems whose solution can be used as a guideline 

^ ■ for its extension. In this respect Supersymmetry (SUSY) [1, 2], and in particular 

O ■ the Minimal Supersymmetric extension of the Standard Model (MSSM) [3-5], is one 

^ ! of the most promising scenarios for physics beyond the Standard Model. In SUSY 

I models the breaking of electroweak symmetry is obtained radiatively at a scale which 

^ I is stabilized by Supersymmetry itself. Moreover, the comparison of MSSM predictions 

and electroweak precision observables provides an overall fit of data [6, 7] which is at 
least as good as that obtained using the SM, and even better in the case of specific 
observables such as (7 — 2 of the muon [8,9]. 

These features render SUSY an appealing framework, and they can explain the 
big effort in the hunting for SUSY not only in the past but also in future years. In 
particular, this search is one of the major goals of the LHC. Experimental studies have 
shown the possibility of early discovery of SUSY with 1 fb~^ of integrated luminosity 
in the inclusive multijet plus missing Et channel [10, 11], provided that the masses 
of squarks and gluinos are not too heavy (e.g. < 2 TeV). Interestingly, the 95 % 
confidence level area of the (jtj-o, "^1/2) plane of the Constrained MSSM (CMSSM) lies 
largely in the region that will be investigated with 1 fb~ at 14 TeV [7]. 
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At hadron colliders, colored particles can be searched for most efficiently. This 
implies that Supersymmetry could be discovered looking at the production of squarks 
and gluinos. Among the others, the process of gluino pair production, 

PP^-g-gX, (1) 

is one of the most important processes leading to the production of colored SUSY 
particles. Indeed, its cross section is large, (9(10 pb) if the gluino mass is (9(600 GeV). 
Moreover, the gluino plays a key role in characterizing SUSY models. The measure 
of the spin of a (supposed to be) gluino [12] and the confirmation of its Majorana 
nature [13-15] would allow not only to distinguish among different beyond Standard 
Model scenarios, but also among MSSM and others SUSY models involving Dirac 
gauginos, such as the N — 1/N — 2 hybrid scheme [16]. 

The total cross section for gluino pair production was computed at Born level 
long time ago [17-20]. NLO SUSY-QCD contributions were computed in Ref. [21]. 
These corrections are positive and large (from 5 to 90%, depending on the masses of 
the squarks and of the gluino), and they reduce appreciably the factorization scale 
dependence. They are included into the pubficly available code Prospino [22]. More 
recent is the resummation of the QCD Sudakov logarithms at the next-to-lcading- 
logarithmic (NLL) accuracy [23,24], and the resummation of the leading Coulomb 
corrections [24]. Their inclusion further stabilizes the prediction against scale variation. 
The NLL contributions are of the order of 2—8% of the NLO QCD predictions, provided 
the squark and gluino masses are 0{1 TeV). In this mass range the contribution of 
the Coulomb corrections amounts up to 5%. The computation of NLO electroweak 
(EW) corrections, of (9(a^a), to the process of hadronic production of a gluino pair 
is still missing. In this paper, we fill this gap computing the full 0{a1a) corrections 
to the process (1) in the framework of the real MSSM. Our computation is part of an 
ongoing project aiming to evaluate the tree-level EW and NLO EW contributions to 
the production of colored SUSY particles at the LHC [25-29]. 

The plan of the paper is the following. In section 2 we briefiy summarize the 
0{a1) contributions to the process (1). In section 3 we describe the partonic processes 
contributing at 0{a1a). Numerical results for the electroweak corrections to gluino 
pair production at the LHC are presented in section 4, while section 5 is devoted to a 
brief discussion on the numerical value of the electroweak corrections at the Tevatron. 
Section 6 summarizes our results. A list of Feynman diagrams is collected in the 
appendix. 

2 Gluino pair production in lowest order 

The leading order contributions to the process (1) are of QCD origin, of 0[a1). At 
lowest order in the perturbative expansion the differential cross section can be written 
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as follows, 



(2) 

with the help of the parton luminosities, defined according to 

fi\p{x)fj\P (^) + fj\p{x)fi\P . (3) 



dLij 1 dx 

dr 1 + dij X . 



fi\H{x) is the momentum distribution of the parton i inside the hadron H. The sum 
runs over the quarks q = u, d, c, s, b. The lower limit on the integral over r. 
To — Amy S, is related to the threshold for the production of the gluino pair. We use 

the convention dax^ to denote the cross section for a partonic process X at a given 
order 0(Q;^a;^) in the strong and electroweak coupling constants. Therefore, da^f_^^^ 

2 

and d(Tgg_^~~ are the lowest order differential cross sections for the partonic processes 

fe) ^ ^(^i) ^(^2), (4) 
9{Pi) 9ip2) g{ki) g{k2), (5) 

respectively. The cross sections are averaged (summed) over the spins and the colors of 
the incoming (outgoing) particles. In this analysis we will consider the five light quarks 
as massless and we approximate the CKM matrix by the unity matrix. We retain the 
mass of the bottom in the Yukawa couplings, owing to the possible enhancement due 
to tan p. We perform our computation in Feynman gauge. 

In lowest order, the partonic cross sections for the processes (4) and (5) can be 
obtained from the Feynman diagrams in Fig. 1 of appendix A. The cross sections can 
be written as 

^^2,0 dt I « >1,0 |2 



where M.]f is the tree-level contribution to the amphtude of the process X. The 
squared amplitudes averaged (summed) over the spins and the colors of the initial 
(final) particles read [17-21] 

El-^«-»r = |"?''1P(2"r' + *^4)+4m|»(^ + ^) 

36(m?s + t1) f 1 1 \ , _ 2 1 , 1 



^99 \ 



A factor 1/2 has been taken into account because of the identical particles in the final 
states. The Mandelstam variables are defined as 



s = (Pl+P2)^ 

t = iP2-hf, 



(8) 
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q,ai 
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q,a — U Tn^g^, 



and u — 2m| — s — t. 



3 0((y.^cx) corrections to the hadronic process 

Gluinos do not interact weakly, thus a pair of ghiinos is neither produced at 0{a'^), from 
qq initial states, nor at 0{asa) via photon-induced processes. Therefore, in contrast 
to squark-anti-squark [25,27,28] and squark-gluino [29] production, EW contributions 
enter only at NLO, and they are at least of 0{ala). The NLO EW contributions to 
the hadronic differential cross section reads as follows. 



(9) 



The g7 and 7g luminosities entering (9) are built according to Eq. (3). Besides the 
virtual corrections and the real photon radiation processes at we consider 

the photon-induced processes leading to the production of a gluino pair together with 
an (anti-)quark. Diagrams and amplitudes are generated with FeynArts [30,31]. The 
reduction of the one-loop integrals is performed with the help of FormCalc [32,33], while 
the scalar one-loop integrals are numerically evaluated using LoopTools. Infrared (IR) 
and coUinear singularities are treated using mass regularization, i.e. giving a small 
mass to the photon and to the five light quarks. 

3.1 qq annihilation with elect roweak loops 

The first class of corrections entering Eq. (9) are the electroweak one-loop corrections 
to the processes (4), yielding the following partonic cross section, 

'''^M = ^' {M^,,Mil.,,} . (10) 

A4q^_,gg is the one- loop electroweak contribution to the amplitude of the quark-anti- 
quark annihilation process. The diagrams responsible for this contribution are dis- 
played in Fig. 4 of appendix A. 
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We treat UV divergences using dimensional reduction. In order to cure the UV 
divergences we have to renormahze the quark and the squark sector at 0{a). Renor- 
mahzation of mass and wavefunction of the quarks and squarks belonging to the first 
two generations, of tan and of the mass of the W boson, has been performed according 
to the procedure described in Ref. [28]. The renormalization of mass and wavefunc- 
tion of the bottom and top quarks and squarks has been widely studied [34-38], and 
several renormalization schemes have been proposed. Each of these schemes has its 
own virtues and drawbacks, we perform our computation using two different renor- 
malization schemes. The first (second) scheme, referred in the following as Rsl (Rs2), 
is the "nib OS" {"nib DR") scheme defined in Ref. [38]. The renormalization of the 
stop-sbottom sector at 0{a) within the Rs2 scheme requires the renormalization of 
the supersymmetric Higgs parameter /j,. This parameters has been defined in the DR 
scheme. 

In the case of bb — >• gg we keep the mass of the b-quark that appears in the couplings. 
In this case, the last twelve diagrams in Fig. 4 of the appendix A have to be considered 
as well. It is well known [39-44] that, in the large tan/3 regime, the tree- level relation 

between the bottom mass nih and the bottom Yukawa couplings yh receives radiative 
corrections that can be strongly enhanced and have to be resummed. Power counting 
in as ta.n (3 shows that the leading tan/5 enhanced contributions, of C(q;^' tan" /3), can 
be accounted for by means of the substitution 



mi 



ml 



mu = 



in the relation between rub and yb- mf'' is the bottom mass in a given renormalization 
scheme, Rs. is defined as 



2a., 



rrig fi tan f3 
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Concerning the Higgs sector, the b — b — h^ coupling is dynamically generated at 0{as). 
This coupling can be enhanced if tan /3 is large and it is worth to include such effects 
modifying the b — 6— Higgs Yukawa couplings. In particular, the effective Lagrangian 
that correctly takes into account these dynamically generated extra-couplings is 



'^Higgs 



m. 



V 

COS a 
cos P 



tan /5 1 - 



A, 



1 + 



tan2/3 
Abtana\ 



tan (3 



sm a 
cos P 



1 - 



A, 



tan a tan (3 



h%b 



H%b 



(13) 
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3.2 qq annihilation with real photon emission 

IR singularities in the virtual corrections are cancelled when the tree-level contribution 
of the partonic process of real photon radiation, 

q{pi) q{P2) ^ g{ki) g{k2) l{h), (14) 

is included. This contribution can be computed using the Feynman diagrams depicted 
in Fig. 2 of appendix A. The integral over the phase space is IR divergent when 
— > 0, while collinear singularities appear whenever • Pi 0. IR and coUinear 
divergences appearing in the phase space integration are regularized using both, phase 
space slicing [45-47] and dipole subtraction [48,49]. The two methods are in good 
numerical agreement, as found also in the case of squark-anti-squark production [28]. 
As already mentioned, IR singularities cancel when the real radiation processes and the 
virtual contributions arc added together, as in Eq. (9). Collinear singularities remain 
and have to be absorbed via the factorization of the parton distribution functions 
(PDFs), c.f. section 3.4. 



3.3 97 and "yq fusion 

The last class of 0{q^^q) contributions to the process (1) are the tree-level contributions 
of the partonic processes 

<l{Pi) l{P2) g{ki) g(k2) qiks), (15) 

7(Pi) fe) g{ki) g(k2) qik^). (16) 

These contributions can be computed from the Feynman diagrams depicted in Fig. 3 
of appendix A. 

Note that, if rriq > nig, the quark in the final state can be the decay product 
of an on-shell squark. If this is the case the last four diagrams depicted in Fig. 3 
become singular. The related poles have to be regularized inserting the width of the 
on-shell squarks into the corresponding propagator. Furthermore, the contribution 
obtained squaring the resonant diagrams has to be subtracted since it arises from the 
production and the subsequent decay of an (anti-)squark through (anti-)quark-photon 
fusion, 

qi M and q ^ g q, 

1 q gq and q ^ g q- (17) 

According to Refs. [21,29], the extraction of the Breit-Wigner pole contribution has 
been performed in the narrow width approximation. 

Collinear singularities arising from initial state emission are again absorbed into 
the PDFs. These singularities are regularized using both, phase space slicing and 
dipole subtraction. The formulae needed can be found in Ref. [45] and in Refs. [49,50], 
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respectively. The results obtained using the two methods agree within the integration 
uncertainty. 

The contribution of this channel is expected to be small owing to the suppression of 
the photon PDF inside the proton. Indeed, the photon PDF is intrinsically suppressed 
with respect to the valence quark PDF by a factor a, since this PDF is originated from 
the emission of a photon from a (anti-)quark. In the SUSY scenarios we consider, the 
contribution of this partonic process amounts up to few percent of the whole O(Qi^a) 
correction. 



3.4 Factorization of initial collinear singularities 

As already mentioned, the universal logarithmic divergences related to the collinear 
splittings 



Q 7, 



Q 



11, 7 



are absorbed into the definition of the PDFs via mass factorization. We factorize the 
(anti-) quark PDFs at 0{a) in the DIS scheme. The effect of this factorization is to 
add the following term into Eq. (9), 



J -.Fact. 
UCTpP^ggX 



2,0 



The functions Ti.'q^ and Hq"' read as follows, 

1 



'TO 

3a 2 



/(2) 



dL 



0.1 I 



dL^ 



(18) 



In 



rrig z{l — z) 



3 1 



21-z 



+ 2z + 3, 



= P,^(^)ln 



tlpl — Z\ 9 

^ - 1 + 8^ - 8^^ 



z 



(19) 



where the splitting functions are 



1 + z' 



Pg^iz) = + (1 - zf. 



Zq is defined as Zq = Amg/irS), while is the charge of the quark q expressed in units 
of the positron charge. The [■■■], distribution is defined as 



C dx [f{x)]+g{x) = f dx f{x)[g{x) - g{l)] - g{l) f dx f{x). 

a J a Jo 



(20) 



In the actual computation, we use the MRST2004qcd parton distribution functions 
at NLO QED and NLO QCD [51]. This fit takes into account QED-effects into 
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parameter 


SPSla' 


SPS2 




250 GeV 


300 GeV 


mo 


70 GeV 


1450 GeV 


Ao 


-300 GeV 





sign(/x) 






tan (5{Mz) 


10.37 


10 



Table 1: MSSM input parameters for the computation of the spectrum of the two 
scenarios considered, mi/2, and Aq are defined at the GUT scale. 



the DGLAP evolution equations and the parametrization of the PDF at the initial 
scale. MRST2004qed PDFs are defined at NLO QCD within the MS mass factoriza- 
tion scheme. As discussed in Ref. [50], the DIS scheme is used for the factorization of 
the 0{a) corrections. 

In our computation we set the renormalization scale, equal to the factorization 
scale, ^F-i and to the gluino mass, i.e. /ir — iif — iTig- 

4 Numerical results, LHC 

For our numerical discussion we use the Standard Model parameters quoted in Ref. [52] . 
The value of the bottom mass in the DR scheme is computed according to Ref. [38] . 
We choose two different SUSY scenarios. The first scenario is the SPSla' suggested by 
the Supersymmetry Parameters Analyses (SPA) [53] project. The second one, called 
SPS2, belongs to the set of Snowmass Points and Slopes, introduced in Ref. [54]. We 
obtain the parameters of the two scenarios with the help of the program SPheno [55], 
starting from the input parameters shown in Table 1. In the SPSla' (SPS2) scenario, 
the gluino mass is 608 GeV (784 GeV). 

4.1 Dependence on the SUSY scenario 

We compute the total hadronic cross section, the results are collected in Table 2. The 
second column shows the lowest order results. The third column shows the sum of the 
lowest order and of the 0{ala) contributions. In the fourth column the contribution 
of the 0{ala) corrections relative to the total result is given, i.e. S is defined as 

EW. NLO 

S = . 

jjQCD, LO _j_ (jEW, NLO 

In the last entry we give an estimate of the statistical error based on an integrated 
luminosity L = 100 fb~^ [56]. We do not distinguish the results in the different renor- 
malization schemes since they agree within the integration error. A priori this is not 
guaranteed. Indeed the Rsl scheme turns out to be unreliable in the scenarios we are 
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point 


^qCD, LO 


^QCD, LO (jEW, NLO 


6 


1 

VL-cr«CD, LO 


SPSla' 


6.1865(6) pb 


6.1822(6) pb 


-0.07% 


0.13% 


SPS2 


1.2127(1) pb 


1.2089(1) pb 


-0.31% 


0.29% 



Tabic 2: Total hadronic cross section for gluino pair production at the LHC (\/S = 
14 TeV). In the second (third) column we show the 0{af) {0{a'j, +a^a)) contribution 
for the points SPSla' and SPS2. In the fourth cohimn the electroweak corrections 
relative to the LO + NLO EW result are given. The last column shows the statistical 
error for an integrated luminosity L = 100 fb~^. 



considering. In this scheme, the finite part of the renormalization constant of the trilin- 
ear coupling, is comparable with the value of the trilinear coupling Ab itself, i.e. 

6A^''/Ab ~ 1, and the perturbative expansion is spoiled. However, the difference among 
the renormalization schemes is as small as few percent of the tree-level bb annihilation 
channel cross section. The latter contributions amount up to several fb, therefore the 
variation of the results in Table 2 is within the integration error. 

As one can see, in the case of the point SPSla' the electroweak corrections are much 
smaller than the statistical uncertainty and so they are not relevant. In the case of the 
point SPS2, the C(q;^q;) corrections are of the same order of the statistical error but 
they are smaller than the theoretical systematic uncertainties such as the uncertainty 
on the PDF parametrization (< 10%) [21] and the factorization scale dependence (from 
3 to 5% if mg < 1 TeV) [24]. 

The invariant mass distribution of the two gluinos is shown in Fig. 5. The EW 
corrections are small, their absolute value being at most of the order of 0.4% of the total 
contribution. Moreover, these corrections do not distort the shape of the distribution. 

In Fig. 6 we consider the distribution of the largest transverse momentum of the 
two gluinos, for brevity we will refer to this observable as "transverse momentum 
distribution". The 0{a'^a) corrections are rather small, the absolute value of their 
contribution relative to the total result is at most 1%, reaching this value in the high 
Pt region, for pt > 1500 GeV. 

4.2 Dependence on the MSSM parameters 

In this subsection we investigate the size of the 0{ala) corrections to gluino pair 
production in a more systematic way, performing a scan over the parameter space 
of the MSSM. The parameters involved in the scan are the independent parameters 
in the second renormalization scheme. We suppose that all the sfermionic soft mass 
parameters are equal and we indicate their value in the Rs2 scheme as Ms^^y. The 
physical masses of the sfermions can be obtained from Mg^gy diagonalizing the mass 
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matrices. Moreover, we consider the surfaces of the parameter space characterized by 
At = At-. With these assumptions there are eight independent parameters involved in 
the scan, namely, 

-^susy, rrig, 11, M2, At, Ab, tan/?, M^o. 

The subregions of the parameter space are chosen imposing the exclusion limits arising 
from SUSY searches at LEP [57] and at the Tevatron [58], and the bound on the mass 
of the light Higgs boson. The physical mass of the light Higgs boson has been computed 
using FeynHiggs 2.5.1 [59-61]. Moreover, each point in the selected regions fulfills the 
condition |Ap[ < 0.025, Ap being the dominant SUSY corrections to the electroweak 
p parameter, corrections arising from top and bottom squarks contributions. 

We perform four different scans. In each scan we select two parameters and we 
study the dependence of the quantity A, 

EW. XLO 

A = 100. 

0-QCD, LO _j_ (jEW, NLO 

We repeat each of these scans for different values of another pair of parameters, while 
the remaining four are fixed to their SPSla' values. Here there is a brief discussion on 
the results of these scans. 

Scan over At and A5 

The results of this scan are displayed in Fig. 7. As expected, A is quite independent 
on the parameter At which enters in the virtual correction of the process bb — > gg and 
in the definition of the mass of the top squarks. This feature is more evident for large 
tan /3 values. A varies only by an amount of the order of few percent for a variation 
of Ab and At over a quite broad range (from —1500 to 1500 GeV ). Note that in the 
whole subregion considered the absolute value of A is of 0(0. 1). 

Scan over tan /3 and Mj^o 

As can be inferred from Fig. 8, the dependence of A on (tan/?, M^o) strongly varies 
for different values of rrig and Msusy As a general result the overall dependence is mild 
for each value of {rrig, Mgu^y)- In all cases the value of |A[ is at most of the order of 0.7. 

Scan over /j, and M2 

As displayed in Fig. 9, A is almost independent on fj, for each value of the pair 
{rrig, Msusy) while the dependence on M2 is more important and particularly pronounced 
when rrig = 1250 GeV and Mg^.y = 730 GeV. In the case of the last three plots the value 
of A is of order —0.3 to —0.1, while in the first plot, characterized by rrig ~ 2 • M^usy, 
the value of A is enhanced for small values of M2 reaching the value of —0.65. 
Notice that the mass of the lightest neutralino and chargino is almost independent on 
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the value of jj, but varies strongly as M2 varies, growing as the value of this parameter 
grows. So this enhancement occurs when charginos/neutralinos are much lighter than 
the gluino. 

Scan over Mgnsy and rUg 

In this scan we investigate the dependence of A on rrig and Mg^^y, which is expected 
to be the most important because of the dependence of the lowest order cross section 
on these parameters. We consider the variation of A as a function of {nig, M^^^y) for 
different values of M2 and tan/9, see Fig. 10. Note that we plot ^ = — A instead of A. 
As a general feature ^ increases as rrig increases and as Mgu^y decreases. The behaviour 
of ^ as a function of Mgu^y and nig is affected by the value of M2 being enhanced for 
smaller values of this parameter. In particular ^ ~ 3 in the region rrig > 1600 GeV, 
Ms^^y < 500 GeV. 

The enhancement of the EW corrections is related to the increasing importance of the 
qq annihilation channel when the production threshold becomes higher. Indeed, the 
minimal value of the parton's momentum fraction rises as the gluino mass rises. Since 
the relative importance of the (anti-)quark PDF increases as the momentum fraction 
of the (anti-)quark increases, the EW corrections grow as the mass of the gluino grows. 
The relative importance of the EW contributions is more pronounced when Ms^sy is 
small owing to the presence of tree-level diagrams with squarks exchanged in the t and 
u channel, c.f. Fig. 1, which are enhanced when the squark masses decrease. 



5 Numerical results, Tevatron 

The EW contributions to gluino pair production are expected to be more important 
at the Tevatron than at the LHC, owing to the enhancement of the quark-ant i- quark 
annihilation channels with respect to the gluon fusion channel. Therefore, it is worth to 
estimate the impact of the EW contributions to gluino pair production at the Tevatron, 
i.e. to the process 

PP^~g~gX. (21) 

The previous analysis can be easily extended to (21), provided that the definition of 
the luminosity, Eq. (3), is replaced by 

For numerical evaluation, we focus on two different points of the MSSM parameter 
space, referred to as TPl and TP2 respectively. These points belong to the region of 
the parameter space of the MSSM used in the data analysis made by CDF and DO 
collaborations [62-64]. We obtain the parameters in these scenarios with the help of 
SPheno, starting from the input parameters at the GUT scale described in Table 3. 



11 



parameter 


TPl 


TP2 


mi/2 


200 GeV 


500 GeV 


mo 


130 GeV 


120 GeV 


Ao 








sign(;u) 






tan /3{Mz) 


3 


3 



Table 3: MSSM input parameters for the computation of the spectrum of the scenarios 
TPl and TP2. 



point 






5 


1 




TPl 


0.16714(1) fb 


0.16691(1) ih 


-0.14% 


61 % 


TP2 


0.048864(3) pb 


0.048256(4) pb 


-1.26% 


3.6 % 



Table 4: Same as Table 2, but considering gluino pair production at the Tevatron, i.e. 
the process PP ggX at VS = 1.96 TeV, and different SUSY scenarios. In this case, 
the last column shows the statistical error for an integrated luminosity L = 2 x 8 fb~^. 



These points are compatible with the experimental limits set by the analysis made by 
the DO collaboration [64] . In particular, the first one corresponds to a scenario in which 
the gluino is heavier than the squarks {rrig ~ 500 GeV, m^_^f ~ 460 GeV), while the 
second one describes a scenario characterized by a light gluino {rrig ~ 340 GeV, and 
m^-_,i ~ 550 GeV). 

In Tabic 4 we show the total hadronic cross section in the two points considered. 
We use the same notation as in Table 2. In the case of the point TPl, the size of the 
electroweak corrections is so small that they will not be visible at the expected final 
integrated luminosity L = 2 x 8 fb~^ . In the case of the point TP2 we obtain a relative 
statistical error of order 4% which is three times bigger than the size of the electroweak 
contributions. Moreover it is worth to notice that the systematic uncertainties affecting 
SUSY searches at the Tevatron are typically greater than 1%. For instance, Ref. [64] 
claims that the jj-p dependence of the total cross section gives an error from 15 to 20%. 

The invariant mass distribution for the two points is shown in Fig. 11. In both 
cases the 0{ala) corrections are small compared to the lowest order results and do 

not change the shape of the distribution. In particular, in the case of the point TPl 
(TP2) EW corrections relative to the total contribution are of the order of —2 to 5% 
(-3 to -1%). 

Similar considerations hold in the case of the transverse momentum distribution, 
shown in Fig. 12. The shape of the distribution is not affected by the insertion of 
the electroweak corrections in both points. The electroweak corrections relative to the 
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total contribution are of order of —2 to 4.5% in the case of the TPl point and of the 
order of —2.5 to —1% in the TP2 scenario. 

6 Conclusions 

In this paper, we have computed the full 0{ala) corrections to gluino pair production 
at the LHC and at the Tevatron. Two different renormalization schemes were used. 
The numerical value of the 0{ala) contribution is rather independent on the renor- 
malization scheme. The treatment of the IR and coUinear singularities was performed 
within two different methods. 

We have studied the numerical impact of the 0{ala) contributions at the LHC in 
two different scenarios and we have performed scans over many regions of the parame- 
ter space. The EW corrections are negative and can be safely neglected. Compared to 
squark-anti-squark [25,28] and squark-gluino [29] production, the EW contributions to 
gluino pair production are less important. The main reason is that the EW contribu- 
tions do not enter the gluon fusion channel, which is the leading tree-level production 
channel in a wide part of the region of the parameter space investigated in this paper. 

We have also provided numerical results for gluino pair production at the Tevatron 
selecting two scenarios belonging to the region of the parameter space investigated 
by the DO and CDF collaborations. Again, the 0{ala) contributions are small and 
negligible. 
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Appendix 



A Feynman Diagrams 

In this appendix wc collect the relevant Feynman diagrams. In the following the label 
S {S^) is used to denote charged (neutral) Higgs bosons. Moreover y° = 7, Z. 




g q 9 



Figure 1: Tree-level diagrams for the processes qq gg and gg gg. 




Figure 2: Tree- level diagrams for the real photon emission process qq — > gg^. 




Figure 3: Trcc-lcvcl diagrams for the process 57 — > ggq. The diagrams for the process 
7? ~^ 990. can be obtained inverting the arrows. 
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Figure 4: One-loop EW diagrams for tlie process qq gg. Diagrams with crossed final 
states are not shown. 



15 



References 



[1] J. Wess and B. Zumino, Supergauge Transformations in Four-Dimensions, Nucl. 
Phys. B70 (1974) 39-50. 

[2] D. V. Volkov and V. P. Akulov, Is the Neutrino a Goldstone Particle?, Phys. 
Lett. B46 (1973) 109-110. 

[3] H. P. Nilles, Super symmetry, Supergravity and Particle Physics, Phys. Rept. 110 
(1984) 1-162. 

[4] H. E. Haber and G. L. Kane, The Search for Super symmetry: Probing Physics 
Beyond the Standard Model, Phys. Rept. 117 (1985) 75-263. 

[5] R. Barbieri, Looking Beyond the Standard Model: The Supersymmetric Option, 
Riv. Nuovo Cim. 11N4 (1988) 1-45. 

[6] J. R. Ellis, S. Heinemeyer, K. A. Olive, A. M. Weber, and G. Weiglein, The 
Supersymmetric Parameter Space in Light of B~ physics Observables and 
Electroweak Precision Data, JHEP 08 (2007) 083, [arXiv: 0706. 0652]. 

[7] O. Buclimueller et. ai. Prediction for the Lightest Higgs Boson Mass in the 
CMSSM using Indirect Experimental Constraints, Phys. Lett. B657 (2007) 
87-94, [arXiv: 0707. 3447]. 

[8] Muon g-2 Collaboration, G. W. Bennett et. ai, Measurement of the Positive 
Muon Anomalous Magnetic Moment to 0.7 ppm, Phys. Rev. Lett. 89 (2002) 
101804, [hep-ex/0208001]. 

[9] Muon g-2 Collaboration, G. W. Bennett et. ai. Measurement of the negative 
muon anomalous magnetic moment to 0. 7-ppm, Phys. Rev. Lett. 92 (2004) 
161802, [hep-ex/0401008]. 

[10] The ATLAS Collaboration, G. Aad et. ai. Expected Performance of the 
ATLAS Experiment - Detector, Trigger and Physics, arXiv: 0901. 0512. 

[11] The CMS Collaboration, Cms physics technical design report, volume ii: 
Physics performance. Journal of Physics G: Nuclear and Particle Physics 34 
(2007), no. 6 995-1579. 

[12] A. Alves, O. Eboli, and T. Plehn, It's a gluino, Phys. Rev. D74 (2006) 095010, 
[hep-ph/0605067]. 

[13] R. M. Barnett, J. F. Gunion, and H. E. Haber, Discovering supersymmetry with 
like sign dileptons, Phys. Lett. B315 (1993) 349-354, [hep-ph/9306204]. 

[14] S. Kraml and A. R. Raklev, Same-sign top quarks as signature of light stops at 
the LHC, Phys. Rev. D73 (2006) 075002, [hep-ph/05 12284]. 



16 



[15] S. Y. Choi, M. Drees, A. Freitas, and P. M. Zerwas, Testing the Majorana Nature 
of Glmnos and Neutralmos, Phys. Rev. D78 (2008) 095007, [arXiv: 0808. 2410]. 

[16] M. M. Nojiri et. ai, Physics Beyond the Standard Model: Supersymmetry, 
arXiv: 0802. 3672. 

[17] G. L. Kane and J. P. Leveille, Experimental Constraints on Gluino Masses and 
Supersymmetric Theories., Phys. Lett. B112 (1982) 227. 

[18] P. R. Harrison and C. H. Llewellyn Smith, Hadroproduction of Supersymmetric 
Particles, Nucl. Phys. B213 (1983) 223. 

[19] E. Reya and D. P. Roy, Supersymmetric particle production at p p collider 
energies, Phys. Rev. D32 (1985) 645. 

[20] S. Dawson, E. Eichten, and C. Quigg, Search for Supersymmetric Particles in 
Hadron - Hadron Collisions, Phys. Rev. D31 (1985) 1581. 

[21] W. Beenakkcr, R. Hopker, M. Spira, and P. M. Zerwas, Squark and gluino 
production at hadron colliders, Nucl. Phys. B492 (1997) 51-103, 
[hep-ph/9610490]. 

[22] W. Beenakker, R. Hopker, and M. Spira, PROSPINO: A program for the 
PROduction of Supersymmetric Particles In Next-to-leading Order QCD, 
hep-ph/9611232. 

[23] A. Kulesza and L. Motyka, Threshold resummation for s quark- antis quark and 
gluino- pair production at the LHC, Phys. Rev. Lett. 102 (2009) 111802, 
[arXiv: 0807. 2405]. 

[24] A. Kulesza and L. Motyka, Soft gluon resummation for the production of 
gluino-gluino and squark- antis quark pairs at the LHC, arXiv: 0905. 4749. 

[25] W. Hollik, M. Kollar, and M. K. Trenkel, Hadronic production of top-squark pairs 
with electroweak NLO contributions, JHEP 02 (2008) 018, [arXiv: 0712. 0287]. 

[26] M. Beccaria, G. Macorini, L. Panizzi, F. M. Renard, and C. Verzegnassi, 

Stop-antistop and sbottom-antisbottom production at LHC: a one-loop search for 
model parameters dependence. Int. J. Mod. Phys. A23 (2008) 4779-4810, 
[arXiv: 0804. 1252]. 

[27] S. Bornhauser, M. Drees, H. K. Dreiner, and J. S. Kim, Electroweak 

Contributions to Squark Pair Production at the LHC, Phys. Rev. D76 (2007) 
095020, [arXiv: 0709. 2544]. 

[28] W. Hollik and E. Mirabella, Squark anti-squark pair production at the LHC: the 
electroweak contribution, JHEP 12 (2008) 087, [arXiv: 0806. 1433]. 



17 



[29] W. Hollik, E. Mirabella, and M. K. Trenkel, Electroweak contributions to 

squark-glmno production at the LHC, JHEP 02 (2009) 002, [arXiv : 0810 . 1044]. 

[30] T. Hahn, Generating Feynman diagrams and amplitudes with FeynArts 3, 
Comput. Phys. Commun. 140 (2001) 418-431, [hep-ph/0012260]. 

[31] T. Hahn and C. Schappacher, The implementation of the minimal 

supersymmetric standard model in FeynArts and FormCalc, Comput. Phys. 
Commun. 143 (2002) 54-68, [hep-ph/0 105349]. 

[32] T. Hahn and M. Perez- Victoria, Automatized one-loop calculations in four and D 
dimensions, Comput. Phys. Commun. 118 (1999) 153-165, [hep-ph/9807565]. 

[33] T. Hahn and M. Ranch, News from FormCalc and LoopTools, Nucl. Phys. Proc. 
Suppl. 157 (2006) 236-240, [hep-ph/0601248]. 

[34] H. Eberl, M. Kincel, W. Majerotto, and Y. Yamada, One-loop corrections to the 
chargino and neutralino mass matrices in the on-shell scheme, Phys. Rev. D64 

(2001) 115013, [hep-ph/0104109]. 

[35] A. Brignole, G. Degrassi, P. Slavich, and F. Zwirner, On the two-loop sbottom 
corrections to the neutral Higgs boson masses in the MSSM, Nucl. Phys. B643 

(2002) 79-92, [hep-ph/020610l]. 

[36] W. Oiler, H. Eberl, W. Majerotto, and C. Weber, Analysis of the chargino and 
neutralino mass parameters at one-loop level, Eur. Phys. J. C29 (2003) 563-572, 
[hep-ph/0304006]. 

[37] W. Hollik and H. Rzehak, The sfermion mass spectrum of the MSSM at the 
one-loop level, Eur. Phys. J. C32 (2003) 127-133, [hep-ph/0305328] . 

[38] S. Heinemeyer, W. Hollik, H. Rzehak, and G. Weiglein, High-precision 

predictions for the MSSM Higgs sector at 0{ai,as), Eur. Phys. J. C39 (2005) 
465-481, [hep-ph/0411114]. 

[39] T. Banks, Supersymmetry and the Quark Mass Matrix, Nucl. Phys. B303 (1988) 
172. 

[40] L. J. Hall, R. Rattazzi, and U. Sarid, The Top quark mass in supersymmetric 
SO (10) unification, Phys. Rev. D50 (1994) 7048-7065, [hep-ph/9306309] . 

[41] R. Hempfling, Yukawa coupling unification with supersymmetric threshold 
corrections, Phys. Rev. D49 (1994) 6168-6172. 

[42] M. S. Carena, M. Olechowski, S. Pokorski, and C. E. M. Wagner, Electroweak 
symmetry breaking and bottom - top Yukawa unification, Nucl. Phys. B426 
(1994) 269-300, [hep-ph/9402253]. 



18 



[43] M. S. Carena, D. Garcia, U. Nierste, and C. E. M. Wagner, Effective Lagrangian 
for the ibH^ interaction in the MSSM and charged Higgs phenomenology, Nucl. 
Phys. B577 (2000) 88-120, [hep-ph/9912516]. 

[44] H. Ebcrl, K. Hidaka, S. Kraml, W. Majcrotto, and Y. Yamada, Improved SUSY 
QCD corrections to Higgs boson decays into quarks and squarks, Phys. Rev. D62 
(2000) 055006, [hep-ph/9912463]. 

[45] V. N. Baier, V. S. Fadin, and V. A. Khoze, Quasireal electron method in 
high-energy quantum electrodynamics., Nucl. Phys. B65 (1973) 381-396. 

[46] A. Denner, Techniques for calculation of electroweak radiative corrections at the 
one loop level and results for W physics at LEP-200, Fortschr. Phys. 41 (1993) 
307-420, [arXiv: 0709. 1075]. 

[47] W. Hollik, T. Kasprzik, and B. A. Kniehl, Electroweak corrections to W-boson 
hadroproduction at finite transverse momentum, Nucl. Phys. B790 (2008) 
138-159, [arXiv: 0707. 2553]. 

[48] S. Dittmaicr, A general approach to photon radiation off fermions, Nucl. Phys. 
B565 (2000) 69-122, [hep-ph/9904440]. 

[49] S. Dittmaier, A. Kabelschacht, and T. Kasprzik, Polarized QED splittings of 
massive fermions and dipole subtraction for non-collinear-safe observables, Nucl. 
Phys. B800 (2008) 146-189, [arXiv: 0802. 1405]. 

[50] K. P. O. Diener, S. Dittmaier, and W. Hollik, Electroweak higher-order effects 
and theoretical uncertainties in deep-inelastic neutrino scattering, Phys. Rev. 
D72 (2005) 093002, [hep-ph/0509084]. 

[51] A. D. Martin, R. G. Roberts, W. J. Stirling, and R. S. Thorne, Parton 
distributions incorporating QED contributions, Eur. Phys. J. C39 (2005) 
155-161, [hep-ph/0411040]. 

[52] Particle Data Group Collaboration, W. M. Yao et. ai. Review of particle 
physics, J. Phys. G33 (2006) 1-1232. 

[53] J. A. Aguilar-Saavedra et. ai, Supersymm,etry parameter analysis: SPA 

convention and project, Eur. Phys. J. C46 (2006) 43-60, [hep-pli/0511344]. 

[54] B. C. Allanach et. ai. The Snowmass points and slopes: Benchmarks for SUSY 
searches, Eur. Phys. J. C25 (2002) 113-123, [hep-ph/0202233] . 

[55] W. Pored, SPheno, a program for calculating supersymmetric spectra, SUSY 
particle decays and SUSY particle production at e+ e- colliders, Comput. Phys. 
Commun. 153 (2003) 275-315, [hep-ph/0301101]. 



19 



[56] E. Accomando, A. Denncr, and C. Meier, Electroweak corrections to W^f and Z'-f 
production at the LHC, Eur. Phys. J. C47 (2006) 125-146, [hep-ph/0509234]. 

[57] ALEPH Collaboration, S. Schael et. ai, Search for neutral MSSM Higgs bosons 
at LEP, Eur. Phys. J. C47 (2006) 547-587, [hep-ex/0602042]. 

[58] DO Collaboration, V. M. Abazov et. al, Search for squarks and gluinos in events 
with jets and missing transverse energy in pp collisions at ■\/s — 1.96-TeV, Phys. 
Lett. B638 (2006) 119-127, [hep-ex/0604029]. 

[59] S. Heinemeyer, W. HoUik, and G. Weiglein, FeynHiggs: a program for the 
calculation of the masses of the neutral CP-even Higgs bosons in the MSSM, 
Comput. Phys. Commun. 124 (2000) 76-89, [hep-ph/9812320]. 

[60] M. Frank et. al, The Higgs boson masses and mixings of the complex MSSM in 
the Fey nman- diagrammatic approach, JHEP 02 (2007) 047, [hep-ph/0611326]. 

[61] G. Degrassi, S. Heinemeyer, W. Hollik, P. Slavich, and G. Weiglein, Towards 
high-precision predictions for the MSSM Higgs sector, Eur. Phys. J. C28 (2003) 
133-143, [hep-ph/02 12020]. 

[62] DO Collaboration, B. Abbott et. al. Search for squarks and gluinos in events 
containing jets and a large imbalance in transverse energy, Phys. Rev. Lett. 83 
(1999) 4937-4942, [hep-ex/9902013]. 

[63] CDF Collaboration, A. A. Affolder et. al, Search for gluinos and scalar quarks 
in pp collisions at ^/s =1.8 TeV using the missing energy plus multijets 
signature, Phys. Rev. Lett. 88 (2002) 041801, [hep-ex/0106001]. 

[64] DO Collaboration, V. M. Abazov et. al.. Search for squarks and gluinos in events 
with jets and missing transverse energy using 2.1 fb~^ of pp collision data at ■\/s 
= 1.96- TeV, Phys. Lett. B660 (2008) 449-457, [arXiv:0712.3805]. 



20 



SPSla^ 




2000 2500 3000 2000 2500 3000 

Mi„„ [GeV] M,„^ [GeV] 



Figure 5: Invariant mass distribution for gluino pair production at the LHC. In the 
left panels we show the LO (black line) and the LO + NLO EW (red line) distribution. 
The two lines are indistinguishable, owing to the smallness of the EW contributions. 
In the right panels the electroweak correction relative to the total result is shown. 
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Figure 7: A = ^ew, nlo/( ^qcd. lo ^ ^ew, nlo ) . ^qo as a function of A and A for 
different values of tan/3 and M^^^y. The other parameters are fixed to their SPSla' 
values. 



Ms„,,=730 GeV, M_=750 GeV 




Figure 8: A = cr 



EW, NLO ^QCD, LO 



100 as a function of tan jS and M^o for 



different values of nig and Msusy The other parameters are fixed to their SPSla' values. 
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Ms„y= 730 GeV.ra =1250 GeV | Ms„„= 730 GeV, IW =750 GeV 




Figure 9: A = ct^^- ^^°/( (t«°°' lo+^ew, nlo ) . ^qq as a function of fi and M2 for different 
values of rrig and Mg^^y. The other parameters are fixed to their SPSla' values. 



tan|3=15, M^=600 GeV | tan p=15, M^=1100 GeV| 




Figure 10: ^ = - A = -a^^' ^° + a^^' ) • 100 as a function of Ms,,, and 

rrig for different values of M2 and tan (3. The other parameters are fixed to their SPSla' 
values. 
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Figure 11: Invariant mass distribution of the two gluinos produced at the Tevatron via 
the process PP — > ggX. In the left panels we show the LO and the EW+NLO EW 
contribution, while the electroweak corrections relative to the total result are shown in 
the right panels. 
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